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ABSTRACT 

We examine the longitudinal distribution of the stars escaping from a cluster along 
tidal tails. Using both theory and simulations, wc show that, even in the case of a star 
cluster in a circular galactic orbit, when the tide is steady, the distribution exhibits 
maxima at a distance of many tidal radii from the cluster. 

Key words: galaxies: kinematics and dynamics - galaxies: star clusters - methods: 
analytical - methods: iV-body simulations 



1 INTRODUCTION 
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Stars escaping from a star cluster in the gravitational field of 
a galaxy form exte nded tidal tails. These have been observed 
in th e Milky Way dGrillmair et alll 19951; iLehmann fc Scholp 



19971: 



Kharchenko. Scholz. fc LehmannT 119971 : 

Combes] l200d: lOdenkirchen et all l200lj ; 
' Lee et al.ll2003l: iBelokurov et all 120061') 
and have o ften been modelled IjCombes, Leon, fc Mevlan| 
19991: Ijohnston, Sigurdsson. fc Hernquistl 1 19991: 



Leon, Mcylan, 
Rockosi et alj|2002 



Yim fc Led 120021 : 
20041: iLee et al l 



Dehnen et al 
T 



2004 



2004: 



di Matteo, Capuzzo Dolcctta, 



Odenkirchen et al 
fc Miocchil 



Koch et all 
20031: 



Capuzzo Dolcetta. Di Matteo. fc Miocchil 



Lee. Lee, fc SunsT 120061: IChumak fc Rastorguevl 



Choi. Weinberg fc Katl 120071 : iFellhauer et all 



Montuori et all 120071 ). They exhibit significant longitudinal 
clumps, or over- and under-densities, which 
to the influence of gravitational 



structure, i.e 
are usually attributed 



shocks, e.g. encounters with spiral arms or giant molecular 
clouds, or passages through a galactic disk or past a galactic 
bulge, or simply motion in a triaxial potential. Here we 
show that clumps also arise when the tidal field is static, in 
the simplest case of a star cluster in a circular orbit about 
an axisymmetric galactic potential. In the next section we 
describe a theoretical approach to the problem, which we 
verify in iV-body calculations in Section 3. The final section 
is a slightly extended summary with discussion. 
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As in standard theory ( Chandrasckharl 1 19421 ) we describe 
the position of a star in a rotating, accelerated coordinate 
system with origin at the centre of the cluster, and axes 
aligned towards the galactic anticentre, in the direction of 
galactic motion of the cluster, and orthogonal to the plane 
of motion of the cluster, respectively. Then the equations of 
motion are 

x - 2fly + {k 2 - 4Q 2 )x = —U x (1) 
y + 2Ux = -U y (2) 
z + u)*z = —U z , (3) 

where Q is the angular velocity of the cluster around the 
galaxy, k is the epicyclic frequency, ui z is the frequency of 
motions orthogonal to the plane of motion of the cluster, 
and U is the gravitational potential of the cluster stars. 

We are interested in escaping stars, and so we shall 
adopt a point-mass approximation for U. We let U = 
—GM/r, where G is the gravitational constant, M is the 
total mass of the cluster and r — y^x 2 + y 2 + z 2 . There are 
equilibria at the Lagrangian points (±X£,0,0), where 

XL ~ Ufi 2 -k 2 ) ' 
The total effective energy of a star (per unit mass) is 

E= \v 2 + ®, 

where v is the star's speed and the effective potential energy 
is 

,2x 2 2 2x GM 
\ )x + W Z Z ) . 

r 



*=~((* 2 



The Lagrangian points are saddle points of where its 
3 GM 

value is $l = . This is also the escape energy from 

2 XL 

the potential well of the cluster in the galactic field. 
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Figure 1. Comparison of eqs. (J4j), ([5j (solid lines) with the solution of eqs. JTJ,j2} (dashed lines), with the point-mass approximation for 
U, with the same initial conditions starting at ±xl (except for a small outward x-component of velocity). The units and tidal parameters 
are such that M = 1000 M and x L = 12.9 pc. 



Except for high-speed escapers created in few-body en- 
counters, stars escape as a result of two-body encounters and 
then an escaper us ually has an energy only slightly above the 
escape energy $ L (|Havlilll970h . Therefore such escapers pass 
close to one of the Lagrange points at slow speed. Thereafter 
the acceleration due to the cluster rapidly diminishes, and 
the motion of the escaping star can be well approximated 
by eqs.(ITJ-Q with the right sides set equal to zero. The so- 
lution starting at time t — at position (xl, 0, 0) with zero 
velocity is then easily found to be 



4fi 2 



hi 



-XL+\1- 



y = 



2Q, ( 4fi 2 . 
1 — XL\smhzt — Kt) 

hi \ K 



XL COS hit (4) 

(5) 

2 = 0. (6) 

The solution for the stars starting at (— Xl,0, 0) can be 
found as easily. 

After leaving the Lagrange point the solution traces an 
oscillatory path along the tidal tail and comes to rest again 
at the cusp of the orbit, which is reached when 

t = to = 2tt/k. (7) 

At this time the location of the star is y = —yc where 



yc = 



47rfi 



4fi 2 



XL. 



Now let us suppose there is a stream of stars following 
this motion. Because the y— component of velocity is zero 
at y — —yc, there is a peak (actually, an infinite peak) in 
the distribution of the values of y at this location. In the 
evolution of a star cluster, however, it will not be visible 
until the first escapers reach this location, at around time 
tc- Furthermore, as the star cluster loses mass by escape, xl 
decreases, and the location of this peak moves closer to the 
cluster as the cluster dissolves. The existence of this peak, 
and the time when it begins to become apparent, are the 
two main predictions of this theory. 

Now we discuss the extent to which the solution we 
have used is a satisfactory approximation to the motion of 
slow escapers from a star cluster. There are two factors to 
consider. 

(i) The effect of neglecting the field of the cluster is shown 
in Fig[T] What will be important in what follows is the y- 
coordinate of the first cusp after the star leaves the Lagrange 
point, and it is clear that this is well approximated. 



(ii) Our other main approximation is the assumption that 
the star leaves the Lagrange point with a very small veloc- 
ity. This cannot be tested thoroughly without recourse to 
an iV-body simulation, which we postpone to Sec(3] But we 
can develop a feel for the result by changing the initial con- 
ditions, e.g. by adding a small initial value to y, of order xl- 
Then the cusp is displaced in y by the same amount, so that 
the relative displacement of the cusp is 



XL 

yc 



4tt Q,{4tt 2 - hi 2 )' 



which is 1/(1271") in the case of a point-mass potential in 
which Q — k. While this is not large, it shows that the den- 
sity enhancement at y ~ —yc will be smoothed out. Similar- 
results are obtained if we vary other initial conditions ap- 
propriately, though usually the cusp disappears, and may 
be replaced by a small epicyclic loop. This calculation also 
shows that the distance of the density enhancement from 
the star cluster is of order 12-7T tidal radii, again in the case 
of a point-mass galactic potential. 

Now we consider the number of stars between the clus- 
ter and the first density enhancement (near y — —yc) - This 
is occupied by stars which escape in a time interval of order 
tc- If we assume that a constant fraction, /x, of stars escapes 
in ea ch half-mass relaxation time (see, for example, ISpitzerl 
Il987h . the relevant number of escaping stars is 

N esc = nN^, 



where 



0.138iV 1/2 r? /2 



G l/2 m l/2 log( 7 AT) ' 

in which rh is the half-mass radius, N is the number of stars, 
m is the (mean) stellar mass, and 7 is a numerical constant 
of order unity. It follows that 

0.138 V rh 

The point of this result is not so much to furnish a numerical 
estimate, but to point out that it is only weakly dependent 
on N. In this sense a single numerical simulation with large 
N is not more useful for investigating the distribution along 
a tidal tail than a single small simulation. However, relax- 
ation proceeds by the cumulative effect of many small en- 
counters, and the mean change (in the energy of a star) per 
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encounter decreases as N increases. Therefore the assump- 
tion of zero velocity at the moment when the star crosses the 
tidal radius becomes increasingly accurate as N increases, 
and we may expect that the density enhancements become 
increasingly pronounced in this limit. 



3 NUMERICAL VERIFICATION 

The numerical simulations used to verify the existence of the 
predicted substructure were carried out with the collisiona l 
iV-body code NBODY6 by Sverre Aarseth j Aarsethl \m$ ) 
on the 52-processor Sun Fire E15k at the Edinburgh Parallel 
Computing Centre. We used the supercomputer as a work 
farm to compute a set of 64 renditions of our test model 
which were then coadded to gain statistical accuracy. 

The test model is a cluster of 1000 stars of each 1 Mq, 
following a Plummer density profile with a half-mass radius 
of 0.8 pc. The cluster orbits about a point mass galaxy of 
M G = 9.565 • 10 10 M Q at a galactic distance of R G = 8.5 
kpc, which corresponds to a rotational velocity of v rot = 220 
kms -1 . 

In this particular c ase of a point-m ass potential the tidal 
radius, x L , is given by (|Spitzerlll987l) 

/ M \ 1/3 
Xl = Rg {3M^) ' 
which yields an initial value of 12.9 pc for the above given 
values. Therefore the cusps of the orbits of escaping stars 
are initially located at ±yc, where 

yc = 12tvxl m 490 pc, 

but this decreases with the ongoing loss of cluster mass. The 
epicyclic frequency of the cluster's orbit about the galactic 
centre is given by 

K = n = Vrot/Ra = 0.026 Myr -1 . 

Hence the time the first stars need to reach this first cusp, 
tc, is about 240 Myr. Before this time we do not expect any 
overdensity within the growing tidal tails. 

The simulations were conducted until total dissolution, 
which occurred at about 2000 Myr. For several times the po- 
sitions of the stars were first centred on the density centre of 
the remaining cluster stars within the tidal radius. This was 
done for each rendition before they were combined. Then 
the number density of stars along the direction of the tidal 
tails, i.e. the y-axis, was counted in bins of 25 pc. 

In Fig. [2] the time evolution of the number density of 
stars along the tidal tails is shown. Due to the symmetry of 
the leading and trailing tail they have been folded together 
to increase the statistical accuracy. As predicted in Section 
2 the first stars need at least the time tc to reach the first 
cusp of their oscillatory orbit along the tail at yc- At 2tc 
the overdensity has built up completely and has become very 
pronounced. As time goes by, the loss of mass causes a de- 
crease in xl, and by eq.([2} the location of the overdensity 
moves closer to the cluster. At Ate (— 1000 Myr) the clus- 
ter mass is only about 300 Mq, hence xl has decreased to 
about 8.6 pc and so j/l is supposed to be at about 320 pc 
at Ate which agrees fairly well with the numerical results. 
Different simulations lose mass at different rates, and the 
resultant distribution of values of xl (and hence yc) con- 
tributes to the smoothing out of the density maxima as t 



increases. Nevertheless, there are signs of density extrema 
also at about ±2yc- 

In Fig. [3] the evolution of the mean speed outside the 
cluster is shown. The velocity change as described in Sec. [5] 
is fairly well reproduced. In the beginning the mean speed is 
dominated by a few high-velocity escapers. But after build- 
ing up, the overdensities correspond to a minimum in the 
mean speed while maxima correspond to underdensities. 
Here again the second peak is visible. Furthermore the initial 
rise of speed verifies the assumption that most stars leave 
the cluster with small velocities. 

Fig. [4] shows a snapshot of the surface density along 
the tidal tails after 1250 Myr. In both the x-y and y-z plane 
the over- and underdensities are fairly pronounced, including 
(especially in the upper diagram) the additional overdensi- 
ties at ±2yc- 



4 CONCLUSIONS 

We have considered the longitudinal structure of the tidal 
tails generated by stars escaping from a star cluster in a 
circular galactic orbit. Even though the tidal field is steady 
(in a rotating frame in which the cluster centre is at rest), 
the tails exhibits density enhancements. These correspond 
to places where escaping stars slow down in their epicyclic 
motion away from the star cluster. Because stars escape with 
a somewhat limited range of initial positions and velocities 
(when they cross the tidal boundary) the location at which 
this happens is similar for most escapers, and the resulting 
density enhancement is visible in simulations provided that 
sufficient numbers of simulations are coadded to improve the 
sampling statistics. We have conducted such simulations, 
using stars of equal mass. The distance of the main clump 
is around 38 tidal radii, in the case of a galaxy modelled as 
a central point mass. 

Similar clumps have been o bserve d by 
ICapuzzo Dolcetta. Pi Matteo. fc Miocchil (|2005l 'l. but 
those examples were a little ambiguous because (i) the 
galactic potential was triaxial, (ii) the orbit was not 
quite circular, and (iii) the code was collisionless. For 
these reasons the mechanism of escape in their model 
was not clear, and could include an unquantified episodic 
contribution from the time-dependence of the tidal field. 
Nevertheless they showed that the clumps coincide with 
places where the escaping stars slow down. The main 
additional contribution of the present paper is that we 
have identified the mechanism causing those overdensities, 
by placing the problem in a more idealised setting and 
providing a theoretical explanation. 

Whether such density enhancements would be visible 
in the tidal tails of real star clusters is complicated by their 
mass spectrum and, especially for globular clusters, by the 
non-circularity of their orbits. Nevertheless our results show 
that the density enhancements in tidal tails, which are ap- 
parently observed, should not necessarily be taken to in- 
dicate time-dependent tides, or encounters with perturbers 
such as giant molecular clouds, dark matter clumps or spiral 
arms. 
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Figure 2. Time evolution of the number density of stars along the tidal tails. Since the leading and trailing tail are symmetric they 
are folded together. Time is given in Myr. As expected the overdensities need more than to to build up; they are most pronounced at 
about 2 — 3 to The predicted location of yc agrees fairly well with the results of the simulations. At itc (~1000 Myr) the cluster has 
lost about 70% of its mass, and hence yc moves closer to the cluster centre. Also the flux of escaping stars diminishes. At 2000 Myr the 
cluster has dissolved almost completely. 
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high-velocity escapers dominate the outer parts. As time goes by the speed profile clearly shows the signature of the epicyclic motion 
described in Sec. [2] (i.e. the dip at about 300pc at t = 1250). Furthermo re the low speed close to the cluster verifi e s the approximation 
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Figure 4. Contour plot of the surface density in the x-y (above) and z-y plane (below). The snapshot shows the superimposed test 
clusters at t = 1250 Myr where the first and second peaks are visible best. Contours correspond to 50, 100, 150, 175, 200, 225 and 250 
stars per bin, where the bin size was 25 pc X 25 pc. 
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